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ABSTRACT 

We make a direct comparison of the derived dark matter (DM) distributions between hydrody- 
namical simulations of dwarf galaxies assuming a ACDM cosmology and the observed dwarf galaxies 
sample from the THINGS survey in terms of (1) the rotation curve shape and (2) the logarithmic inner 
density slope a of mass density profiles. The simulations, which include the effect of baryonic feedback 
processes, such as gas cooling, star formation, cosmic UV background heating and most importantly 
physically motivated gas outflows driven by supernovae (SNe), form bulgcless galaxies with DM cores. 
We show that the stellar and baryonic mass is similar to that inferred from photometric and kinematic 
methods for galaxies of similar circular velocity. Analyzing the simulations in exactly the same way 
as the observational sample allows us to address directly the so-called "cusp/core" problem in the 
ACDM model. We show that the rotation curves of the simulated dwarf galaxies rise less steeply than 
CDM rotation curves and are consistent with those of the THINGS dwarf galaxies. The mean value 
of the logarithmic inner density slopes a of the simulated galaxies' dark matter density profiles is 
~— 0.4 ± 0.1, which shows good agreement with a=— 0.29 ± 0.07 of the THINGS dwarf galaxies. The 
effect of non-circular motions is not significant enough to affect the results. This confirms that the 
baryonic feedback processes included in the simulations are efficiently able to make the initial cusps 
with a ~— 1.0 to —1.5 predicted by dark-matter-only simulations shallower, and induce DM halos 
with a central mass distribution similar to that observed in nearby dwarf galaxies. 



Subject headings: Galaxies: dark matter halos - galaxies: 
galaxies: kinematics and dynamics 
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1. INTRODUCTION 

The dark matter (DM) distributions at the centers of 
galaxies have been intensively discussed from both ob- 
servational and theoretical sides for almost two decades 
ever since high-resolution N-body dark matter simula- 
tions assuming a universe dominated by cold dark matter 
(CDM) and a cosmological constant A were performed. 
The ACDM simulations have invariably predicted a dark 
matter density distribution which diverges toward the 
centers of galaxies (Moore 1994; Navarro, Frenk & White 
1996, 1997; Moore et al. 1999b; Ghigna et al. 2000; 
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Klypin et al. 2001; Power et al. 2002; Stoehr et al. 2003; 
Navarro et al. 2004; Reed et al. 2005; Diemand et al. 
2008). In order to describe such cusp-like dark matter 
distributions, Navarro, Frenk and White (1995; 1996) 
proposed a profile (hereafter NFW profile) which can be 
approximated by two power laws, p ~r _1 and p r^r~ 30 
to describe the inner and outer regions of a dark matter 
halo, respectively. In particular, the central cusp feature 



with pr- 



-1.0 



has provided a useful test for ACDM cos- 



mology and sparked interest in seeking constraints by ob- 
serving mass distributions at the center of galaxies (e.g., 
McGaugh et al. 2003). It is worth mentioning that re- 
cent simulations show shallower slopes with a ~— 0.8 at 
radii of 120 pc (Stadel et al. 2009) although DM slopes 
at r < 1 kpc are generally steeper (see also Navarro et al. 
2010). 

The prediction of a central cusp from the ACDM model 
has been seriously challenged by observations of dwarf 
and Low Surface Brightness (LSB) disk galaxies. Obser- 
vations of dwarf and LSB galaxies generally indicate a 
constant matter distribution toward their centers, with 
mass density profiles with a kpc sized core radius. This 
discrepancy of central dark matter distributions in dwarf 
galaxies in ACDM simulations and observations is re- 
ferred to as the "cusp/core" problem. This is a funda- 
mental problem for ACDM together with the likely con- 
nected substructure and angular momentum problems 
(Moore et al. 1999a; Klypin et al. 1999; Simon & Geha 
2007; Navarro & White 1994; Dutton 2009). 

Compared to other types of galaxies, dwarf galaxies 
provide us with a good opportunity for measuring the 
dark matter distribution near the centers of galaxies due 
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to the fact that they have a simple dynamical struc- 
ture (disk galaxies without bulges) but also have low 
baryon fractions and hence less of a dynamical contribu- 
tion by baryons (de Blok & McGaugh 1997). Recently, 
Oh et al. (2011) presented high-resolution dark matter 
density profiles of 7 dwarf galaxies taken from "The Hi 
Nearby Galaxy Survey" (THINGS; Walter ct al. 2008). 
The high quality data from THINGS significantly min- 
imizes observational uncertainties and thus allows us to 
investigate the central dark matter distribution of the 
dwarf galaxies in detail. Mass models of stars and gas 
are constructed using the Spitzer IRAC 3.6^m data from 
the "Spitzer Infrared Nearby Galaxies Survey" (SINGS; 
Kennicutt et al. 2003) and the total integrated THINGS 
Hi map, respectively. The kinematics of baryons are then 
subtracted from the total kinematics of the galaxies in 
order to derive the dark matter distributions. 

One of the main results of Oh et al. (2011) was to ro- 
bustly confirm that the rotation curves of the 7 THINGS 
dwarf galaxies rise too slowly to be consistent with 
a cusp feature at their centers. Moreover, the mean 
value of the logarithmic inner slopes of the dark mat- 
ter density profiles is a = —0.29 ± 0.07, which signifi- 
cantly deviates from the a ~— 1.0 predicted from dark- 
matter-only simulations. The exquisite data used in 
this study allowed unprecedented treatment of the ef- 
fects of observational uncertainties, such as beam smear- 
ing, center offset and non-circular motions, which may 
play a role in hiding central cusps (Blais-Oucllcttc et al. 
1999; van den Bosch et al. 2000; Bolatto ct al. 2002; 
Swaters et al. 2003; Simon et al. 2003; Rhee et al. 2004; 
Gentile et al. 2005; Spckkens et al. 2005; Oh et al. 
2008). The results have thus significantly strengthened 
the observational evidence that the dark matter distri- 
bution near the centers of dwarf galaxies follows a near- 
constant density core. 

Before using these results as a repudiation of cold dark 
matter, however, one must also examine carefully the 
modeling on which the central cusp predictions are based, 
which are usually done using N-body simulations that 
include only the effects of gravity on structure forma- 
tion. Although baryons make up only ~14% of the mat- 
ter of the Universe, this dissipative constituent of the 
Universe cools with cosmological time and accumulates 
in the central regions of DM halos, making up a dynami- 
cally important fraction. Several mechanisms have been 
proposed whereby these central baryons can affect the 
central cusp-like dark matter distribution which is found 
in pure dark matter simulations. A rapid change in po- 
tential (faster than the dynamical time) due to star burst 
triggered outflows, is one mechanism which has been 
shown to be capable of transforming cusp-like profiles 
into (flatter) cores (Navarro et al. 1996; Read & Gilmorc 
2005 although see Gnedin & Zhao 2002). Supernova 
driven random bulk motions of gas in proto-galaxies 
(Mashchenko et al. 2006) has also been shown in mod- 
els to flatten cusps, as have the effects of dynamical 
friction acting on gas clumps (El-Zant et al. 2002), and 
the transfer of angular momentum from baryons to the 
dark matter (Tonini et al. 2006). Modeling an inhomo- 
geneous multi-phase, interstellar medium is critical for 
simulating the baryonic feedback processes in galaxies 
(Robertson & Kravtsov 2008; Ceverino & Klypin 2009). 
Yet cosmological simulations have not, until recently, 



been able to achieve enough resolution to model even 
such an inhomogeneous ISM, but have been forced to 
treat the important processes of star formation and feed- 
back as "sub-grid" physics, averaging star formation and 
supernova feedback over large volumes, compared to the 
typical structural scales (<1 kpc) of small galaxy disks. 

Most recently, Governato et al. (2010) have performed 
high-resolution cosmological N-body+Smoothed Parti- 
cle Hydrodynamic (SPH) simulations of dwarf galaxies 
under the ACDM paradigm, that include the effect of 
baryonic feedback processes, such as gas cooling, star 
formation, cosmic UV background heating and most im- 
portantly physically motivated gas outflows driven by 
SNe. The major finding of Governato et al. (2010) was 
that once star formation is associated with high den- 
sity gas regions, a significant amount of baryons with 
low angular momentum is efficiently removed by strong 
SNe-driven injection of thermal energy and the follow- 
ing gas outflows that carry an amount of gas at a rate 
of 2-6 times the local star formation rate. The large 
scale outflows in turn induce two effects: the loss of 
low angular momentum gas from the central regions pre- 
vents the formation of bulges in low mass systems. Sec- 
ondly, the clumpy nature of the gas and rapid ejection 
on short timescales has dynamical effects on the dark 
matter potential, creating a shallower density profile 
(see Mashchenko et al. 2008; Ceverino & Klypin 2009; 
Mo & Mao 2004; Mashchenko et al. 2006). 

Moreover, the simulated galaxies have a z=0 baryonic 
budget consistent with photometric and kinematic esti- 
mates (van den Bosch et al. 2001; McGaugh et al. 2010). 
The kinematic properties of the simulated dwarf galaxies 
are very similar to those of the THINGS dwarf galaxies in 
terms of their maximum rotation velocity (~60 kms -1 ) 
and dynamical mass (~10 9 M ), allowing a direct com- 
parison between the simulations and observations to be 
made. 

A first estimate of the bulge-to-disk ratio and 
DM properties of these models was presented in 
Governato et al. (2010) by fitting a Sersic profile to ar- 
tificial i-band images and fitting a rotation curve to the 
rotational motions of cold gas using a tilted-ring analy- 
sis. Here, we make more direct comparisons with obser- 
vations, applying the same analysis techniques and tools 
to the simulations as done in the most recent observa- 
tional sample of analoguous galaxies. The determina- 
tion of the DM slopes for the simulations has to be done 
in the same way as in observations since there are var- 
ious aspects that can affect the answer, such as beam 
smearing, center offset and non-circular motions, which 
were not taken into account when the slope was derived 
from the raw simulation data in Governato et al. (2010). 
This approach will have two crucial benefits: (1) pro- 
vide a strong test of theoretical predictions and (2) test 
the extent of observational biases possibly associated to 
current models of baryon/DM mass decompositions in 
real galaxies, in particular the recovery of non-circular 
motions and pressure support induced by SN feedback 
(Valenzuela et al. 2007; Dalcanton & Stilp 2010). 

The structure of this paper is as follows. The sim- 
ulations are described in Section 2. In Section 3, we 
present the dark matter mass modeling of the simulated 
dwarf galaxies. Section 4 compares the derived dark mat- 
ter distributions from the simulations with those of the 
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FlG. 1. — Gas and stellar components of DG1. (a): Simulated integrated gas map (moment 0). Contours start from a column density of 
10 20 cm -2 to 10 22 cm - 2 in steps of 10 cm~~ 2 . (b): Hermite hs velocity field. Velocity contours run from —50 km s~ 1 to 50 kms — with a spacing 
of lOkms - , (c): Velocity dispersion map (moment 2). Velocity contours run from 5kms _1 (gray) to 15kms — 1 (black) with a spacing of 4kms — 1 . 
(d)(e)(f):: Total intensity maps in simulated optical 5, R and Spitzer IRAC 3.6 fim bands with superimposed contours of the gas map in the panel 
(a). 

10 10 M Q (the virial mass is measured within the virial 
radius R V i r , the radius enclosing an overdensity of 100 
times the cosmological critical density). 

To evolve the simulations described here we have 
used the fully parallel, N-body+SPH code GASOLINE 
to compute the evolution of both the collisionless and 
dissipativc component in the simulations. A detailed 
description of the code is available in the literature 
(Wadsley et al. 2004). The version of the code used 
in this paper includes radiative cooling and accounts 
for the effect of a uniform background radiation field 
on the ionization and excitation state of the gas. The 
cosmic ultraviolet background is implemented using the 
Haardt-Madau model (Haardt & Madau 1996), includ- 
ing photoionizing and photoheating rates produced by 
Pop III stars, QSOs and galaxies starting at z = 9. 
We use a standard cooling function for a primordial 
mixture of atomic hydrogen and helium at high gas 
temperatures and we include low temperature cooling 
(Mashchenko et al. 2006). 

In the simulations described in this paper star forma- 
tion occurs when cold gas reaches a given threshold den- 
sity (e.g., Stinson et al. 2006) typical of actual star form- 
ing regions (we used 100 atomic-mass-unit (amu) cm -3 ). 
Star formation (SF) then proceeds at a rate proportional 
to pla S , i-e. locally enforcing a Schmidt law. The adopted 



THINGS dwarf galaxies. Lastly, we summarize the main 
results of this paper and conclusions in Section 5. 

2. THE SPH TREECODE GASOLINE 

The halos for these simulations were selected from a set 
of large-scale, low-resolution, dark matter only simula- 
tion run in a concordance, flat, A-dominated cosmology: 
n = 0.24, A=0.76, h = 0.73, ct 8 = 0.77, and n b = 0.042 
(Verde et al. 2003). The size of the box, 25 Mpc, is large 
enough to provide realistic torques for the small galaxies 
used in this work. The power spectra to model the ini- 
tial linear density field were calculated using the CMB- 
FAST code to generate transfer functions. To include 
the effects of cosmic torques from the large scale struc- 
ture we used the volume renormalization (or "zoom in" ) 
technique (Katz & White 1993). Dark matter particle 
masses in the high resolution regions are 1.6 xlO 4 M Q , 
while the mass of star particles is only 1000 Mq and the 
force resolution, i.e., the gravitational softening, is 86 
pc. In total, at z=0 there are 3.3 x 10 6 particles within 
the virial radius of the simulated dwarf galaxy, hereafter 
referred to as DG1. For all particle species, the gravita- 
tional spline softening, e(z), was evolved in a comoving 
manner from the starting redshift (z ~ 100) until z=8, 
and then remained fixed at its final value from z=8 to 
the present. At z = the virial masses of the halos that 
we studied in this paper are 3.5 (DG1) and 2.0 (DG2) x 
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feedback scheme is implemented by releasing thermal en- 
ergy from SNe into the gas surrounding each star particle 
(Stinson et al. 2006). The energy release rate is tied to 
the time of formation of each particle (which effectively 
ages as a single stellar population with a Kroupa IMF). 
To model the effect of feedback at unresolved scales, the 
affected gas has its cooling shut off for a time scale pro- 
portional to the Sedov solution of the blast wave equa- 
tion, which is set by the local density and temperature 
of the gas and the amount of energy involved. The effect 
of feedback is to regulate star formation in the discs of 
massive galaxies and to greatly lower the star formation 
efficiency in galaxies with peak circular velocity in the 50 
< V c <150 kms" 1 range (Brooks et al. 2007). At even 
smaller halo masses (V c < 20-40 kms -1 ) the collapse 
of baryons is largely suppressed by the cosmic UV field 
(Quinn et al. 1996; Okamoto et al. 2008; Gnedin 2010). 

Other than the density threshold only two other pa- 
rameters are needed, the star formation efficiency (eSF 
= 0.1) and the fraction of SN energy coupled to the ISM 
(eSN=0.4). The model galaxies studied in this paper are 
those published in Governato et al. (2010). We verified 
that the addition of full metal cooling and the increase 
of eSN to 1 does not substantially change the structural 
properties of the galaxy. However the amount of stars 
formed decreases by 40%. 

As a benchmark against which the effects of bary- 
onic feedback processes can be gauged, Governato et al. 
(2010) performed additional runs for analoguous model 
galaxies of DG1, called DG1DM and DG1LT. DG1DM 
uses the same initial conditions as DG1 but includes 
only the dark matter component. It has a DM slope of 
a ~— 1.3 similar to those found from similar simulations 
(Springel et al. 2008). DG1LT is a version of DG1 using 
a lower density threshold (0.1 amucm -3 ) at a lower res- 
olution (the force resolution is ~116 pc), and outflows in 
this model are negligible. DG1LT has a cusp and its DM 
slope is similar to that of DG1DM (a ^ — 1.3). 

To properly compare the outputs from the sim- 
ulation to real galaxies and make accurate esti- 
mates of the observable properties of galaxies (e.g., 
Sanchez- Janssen et al. 2010), we used the Monte Garlo 
radiation transfer code SUNRISE (Jonsson et al. 2010) 
to generate artificial optical images and spectral energy 
distributions (SEDs) of the outputs of our run. SUN- 
RISE allows us to measure the dust reprocessed SED of 
every resolution element of the simulated galaxies, from 
the far UV to the far IR, with a full 3-D treatment of 
radiative transfer. We place the simulated galaxies at 
an inclination of 60 degrees. However, in applying our 
analysis tools, the inclination is considered a free param- 
eter, in keeping with the techniques applied to observed 
galaxies. Filters mimicking those of the SDSS survey are 
used to create mock observations. 

We note that in the runs adopting the "high thresh- 
old" SF, feedback produces winds that are comparable 
in strength to those found in real galaxies of similar 
mass. However, in our simulations the cold ISM is still 
only moderately turbulent (~ 5-10 kms -1 at z = 0), 
consistent with observations, and the galaxies match 
the observed stellar and baryonic Tully-Fisher relation 
(Governato et al. 2009), as the SF efficiency is regulated 
to form an amount of stars similar to that of real dwarf 
galaxies of similar rotation velocity. 



In the following sections, we perform dark matter mass 
modeling of both DG1 and DG2 in exactly the same 
way as the THINGS dwarf galaxies sample described in 
Oh et al. (2008, 2011). 

3. THE MASS MODELING OF THE SIMULATED DWARF 
GALAXIES 

3.1. The rotation curves 

We first construct the data cubes of DG1 and DG2 by 
tracing the motions of the gas component. The beam and 
velocity resolutions of the cubes are ^6" (corresponding 
to 100 pc at a distance of 4 Mpc) and 2.0 kms -1 , respec- 
tively. As an example, we show the integrated gas map, 
velocity field and velocity dispersion map extracted from 
the cube for DG1 in Fig. 1. For the velocity field, we use 
the Gauss-Hermite polynomial to model the skewness of 
a non-Gaussian profile caused by multiple velocity com- 
ponents (van der Marel & Franx 1993). This function 
includes an extra parameter, called /13, that measures 
the skewness of the Gaussian function, and thus provides 
more reliable central velocities even for profiles with sig- 
nificant asymmetries. As discussed in Oh et al. (2011) 
(see also dc Blok ct al. 2008), the hermite h 3 velocity 
field gives a robust estimate for the underlying circular 
rotation of a galaxy in which non-circular motions are 
insignificant, like DG1 and DG2. Hermite polynomi- 
als have also been used to extract the velocity fields of 
the THINGS galaxies sample (de Blok ct al. 2008). 

Oh et al. (2011) use the bulk velocity fields when de- 
riving the rotation curves of the THINGS dwarf galaxies, 
except for DDO 154 and M81dwB in which non-circular 
motions are insignificant. Compared to other types of 
velocity fields (e.g., intensity- weighted mean, hermite /13, 
single Gaussian and peak velocity fields) , the bulk veloc- 
ity field more effectively minimizes the effect of small- 
scale random motions on the derived kinematics of a 
galaxy (see Oh et al. 2008 for details). However, for a 
galaxy that is not significantly affected by non-circular 
motions like DG1 and DG2, the bulk velocity field is 
nearly identical to the hermite /13 velocity field. 

In Fig. lb, the iso-velocity contours of the extracted 
hermite /13 velocity field are distorted in some regions, 
suggestive of non-circular motions. These are mainly due 
to the SN-driven gas outflows in the simulations. How- 
ever, despite the presence of non-circular motions the 
overall pattern of the galaxy rotation is well recovered in 
the velocity field. 

It is customary to use a set of concentric tilted-rings 
to model the velocity field of a galaxy, each with its 
own kinematic center (XPOS, YPOS), inclination INCL, 
position angle PA, expansion velocity VEXP, systemic 
velocity VSYS and rotation velocity VROT (Begcman 
1989). PA is the angle measured counter-clockwise from 
the north direction in the sky to the major axis of the 
receding half of the galaxy. 

Assuming an infinitely thin disk, we fit these tilted- 
rings to the hermite /13 velocity fields of DG1 and DG2 
to derive their rotation curves. The derived tilted-ring 
models for DG1 and DG2 are shown in Fig. 2. The er- 
rorbar in the rotation velocities indicates the dispersion 
of individual velocity values found along a tilted-ring. 
As mentioned earlier, gas outflows driven by SN explo- 
sions in the simulations cause non-circular motions at 
some regions and induce larger uncertainties in the fit- 
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FlG. 2. — The tiltcd-ring models derived from the hermite /13 velocity fields of DG1 (upper) and DG2 (lower). The open circles in all panels 
indicate the fit made with all parameters free. The filled black circles and solid lines in all panels show the finally adopted tiltcd-ring models as a 
function of galaxy radius. The dashed lines indicate the geometrical parameters used when extracting the rotation velocities (dotted lines in the 
VROT panels) from the true mass distributions of the simulated galaxies. 



ted tiltcd-rings. These are seen as the scatter of the very 
first run results (open circles) with all ring parameters 
free in Fig. 2. However, the local scatter averages out af- 
ter several iterations and the final rotation curves (solid 
lines) in Fig. 2 seem to give a good description of the 
underlying kinematics of DG1 and DG2. As shown in 
the VROT panels of Fig. 2, this can be confirmed by 
the true rotation velocities (dotted lines) derived using 
the full three dimensional mass distributions rather than 
the projected two dimensional surface density profiles of 
the simulated galaxies. Despite not only the uncertain- 
ties but also the assumption of the tilted-ring analysis, 
i.e., an infinitely thin disk, the difference between the 



true and derived rotation velocities is less than 5 kms -1 , 
equivalent to about twice the velocity resolution of the 
cubes. 

For galaxies whose velocity dispersions are large 
enough compared to their maximum rotation velocities, 
we need to correct for the asymmetric drift to obtain 
more reliable rotation velocities (Bureau & Carignan 
2002). However, the 2 nd moment maps of DG1 and DG2 
show small velocity dispersions (~7 kms -1 ) compared 
to the maximum rotation velocities (~60 kms -1 ), and 
the pressure support is insignificant with respect to the 
circular rotations. Therefore, we ignore the asymmetric 
drift corrections for DG1 and DG2. 
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FlG. 3. — Mass models for the baryons of DG1 (t°p) and DG2 (bottom), (a): Azimuthally averaged surface brightness profiles (corrected for 
inclination) in the simulated 3.6/iin, R, V, and B bands (top to bottom) derived applying the tilted-ring parameters shown in Fig. 2. Note that 
the i?-band surface brightness profile is shifted upward by 0.5 mag for clarity. The dashed line for the 3.6/im profile indicates a least-squares fit 
to the data, the radial range over which the fit is made being indicated by the filled circles, (b): The T* in the 3.6/im band derived from stellar 
population synthesis models. The short and long dashed lines show the T^" 6 values derived using optical colors B — V and B — i?, respectively. The 
solid line indicates the mean value adopted as the final T^' 6 . (c): The stellar mass surface density derived from the 3.6/xm surface brightness in (a) 
using the T^" 6 value shown in the panel (b). The red open circles indicate the true profile derived from the simulations, (d): The rotation velocity 
for the stellar component derived from the stellar mass density profile (dots) in the panel (c). The dashed line shows the true rotation velocity for 
the stellar component derived from the true profile (open circles) shown in the panel (c). (e): The radial mass surface density distribution of the 
gas component scaled by 1.4 to account for He and metals. The red open circles indicate the true profile derived from the simulations, (f): The 
gas rotation velocity derived from the gas surface density profile (dots) in the panel (c). The dashed line shows the true rotation velocity for the 
gas component derived from the true profile (open circles) shown in the panel (e). 
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TABLE 1 

Properties of the simulated and the THINGS dwarf galaxies 



Name 


D 


(Incl.) 


20 


M B 


Vmax 


Rmax 


M dyn 


Mhalo 


M s tar 


M gas 




(Mpc) 


n 


(kpc) 


(mag) 


(km S — 1 ) 


(kpc) 


(1O 9 M ) 


(10 9 M Q ) 


(10 8 M ) 


(10 s M @ ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


IC 2574 


4.0 


55.7 


0.57 


-18.1 


77.6 


10.4 


14.6 


53.2 


10.38 


18.63 


NGC 2366 


3.4 


39.8 


0.34 


-17.2 


57.5 


5.6 


4.3 


76.9 


2.58 


6.98 


Holmbcrg I 


3.8 


13.9 


0.55 


-14.8 


38.0 


1.5 


0.5 


33.1 


1.25 


2.06 


Holmbcrg II 


3.4 


49.6 


0.28 


-16.9 


35.5 


7.1 


2.1 


4.3 


2.00 


7.41 


M81 dwB 


5.3 


44.8 


0.09 


-14.2 


39.8 


0.8 


0.3 


870.9 


0.30 


0.31 


DDO 53 


3.6 


27.0 


0.14 


-13.4 


32.4 


2.0 


0.5 


2.1 


0.18 


0.85 


DDO 154 


1.3 


66.0 


0.20 


-14.2 


53.2 


8.2 


5.1 


2.4 


0.26 


3.58 


DG1 


4.0 


60.0 


0.41 


-15.9 


60.0 


5.0 


4.2 


35.0 


1.81 


8.04 


DG2 


1.0 


60.0 


0.28 


-15.6 


56.4 


5.0 


3.7 


20.0 


0.80 


5.92 



Note. — Distance as given in Walter et al. (2008). For DG1 and DG2, we assume they are at a distance of 4 

Mpc; (2): Average value of the inclination derived from the tiltcd-ring analysis; (3): The vertical scale height of disk; (4): 
Absolute B magnitude as given in Walter et al. (2008); (5): Maximum rotation velocity; (6): The radius where the rotation 
velocity V ma x of the flat part of the rotation curve is measured; (7): Dynamical mass from the measured V ma x and R max - 
(8): Halo mass A/200 determined from V200 using Eq. 3. For the THINGS dwarf galaxies we use the V200 values fitted using 
c fixed to 5. See Section 4.1 for more details. For DG1 and DG2, we show the virial mass measured within the virial radius 
R v i r enclosing an overdensity of 100 times the cosmological critical density; (9): Stellar mass derived in Section 3.2. The 
stellar mass of DDO 154 is from de Blok et al. (2008); (10): Gas mass derived in Section 3.2. The gas mass of DDO 154 is 
from Walter ct al. (2008). 
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FlG. 4. — The Disk-halo decompositions of the rotation curves of DGl (left) and DG2 (right). The gray dots represent the total rotation 
curves derived from the hcrmite /13 velocity fields. The dotted and dash-dotted lines indicate the rotation curves of the gas and stellar components, 
respectively. The open circles show the dark matter rotation curves derived subtracting the rotation curves of baryons from the total rotation curves. 
The dashed and solid lines show the best fitted NFW and pseudo-isothermal halo models to the dark matter rotation curves, respectively. The 
reduced x 2 value for each halo model is denoted on the panels. The lower panels show the velocity residuals between the dark matter rotation curves 
and the best fitted halo models. The dots and open circles indicate the results from the NFW and pseudo-isothermal halo models, respectively. 
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3.2. The mass models of baryons 

The derived rotation curves in the previous section rep- 
resent the total kinematics of the galaxies, including not 
only dark matter but also the stellar and gas components. 
We therefore construct mass models of the baryons and 
subtract them from the total kinematics to separate the 
dark matter component only. 

For the gas component, we first derive the gas surface 
density profile by applying the derived tiltcd-rings in Sec- 
tion 3.1 to the integrated gas map shown in the upper- left 
panel of Fig. 1. The derived gas surface density profile is 
then scaled by a factor of 1.4 to account for helium and 
metals. The resulting gas surface density profiles of DG1 
and DG2 are given in panels (e) of Fig. 3. From these, 
we calculate the rotation velocity due to the gas compo- 
nent, assuming an infinitely thin disk. The derived gas 
rotation velocities of DG1 and DG2 arc shown in panels 
(f) of Fig. 3. We also overplot the true values derived 
using the full three dimensional mass distribution of the 
gas components of DG1 and DG2 as shown in the open 
circles in the panels (e) of Fig. 3. For DG1, the true and 
derived gas surface density profiles are similar but the 
derived rotation velocity is systematically higher than 
the true one. This can be due to the assumption of "an 
infinitely thin disk 1 " which makes one overestimate the 
gas rotation velocity of a galaxy with a considerable gas 
thickness. However, the velocity difference (~5 kms -1 ) 
is not very significant. For DG2, the true gas velocity is 
systematically higher in the range 0-3 kpc, particularly 
in the inner region. Likewise, the true surface density 
profile is higher than the derived one in the range 0- 
3 kpc. As shown in Fig. 2, this is partially because the 
smaller inclination value (60°) is used for extracting the 
true gas surface density profile from the simulation in the 
inner region of DG2. In addition this may also be due to 
significant vertical gas outflows or flaring driven by SNc 
perpendicular to the disk of the galaxy. 

Similar to the gas component, we derive the surface 
density profiles for the stellar components of DG1 and 
DG2. For this, as shown in Fig. 1, we use the simulated 
B, V , R and Spitzer IRAC1 3.6/xm images. In particu- 
lar, the Spitzer IRAC1 3.6/xm image is useful for tracing 
the underlying old stars which are usually the dominant 
stellar population in dwarf galaxies. For the same rea- 
son, Spitzer IRAC 3.6/im images have also been used for 
making mass models for the stellar components of the 
THINGS galaxies sample (de Blok et al. 2008; Oh et al. 
2008, 2011). 

We derive the surface brightness profiles of DG1 and 
DG2 by applying the tilted-rings derived in Section 3.1 
to their simulated B, V and R as well as 3.6/xm images. 
Both DG1 and DG2 are bulgeless as found from fitting 
a Scrsic profile to their i band images (Governato et al. 
2010). To convert the surface brightness profiles to the 
mass density profiles in units of M. & pc~ 2 , we obtain 
the 3.6/zm mass-to-light (T^ 6 ) values using an empir- 
ical relation between and optical colors based on 
the Bruzual & Chariot (2003) stellar population synthe- 
sis models (Oh et al. 2008; see also Bell & de Jong 2001). 

1 The gas rotation velocities assuming "an exponential density 
law" with scale heights in the range of 0.5-2.5 kpc are somewhat 
similar to the true one but slightly higher in the outer regions, 
possibly due to SNc-drivcn gas outflows or flaring. 



The T^- 6 values used for DG1 and DG2 are shown as the 
solid lines in panels (b) of Fig. 3, and the resulting stellar 
surface density profiles are given in panels (c) of Fig. 3. 
The true surface density profiles for the stellar compo- 
nents of DG1 and DG2 are also overplotted as the open 
circles in Fig. 3, and they are similar to the derived ones. 
This can be treated as circumstantial evidence that the 
assumption used for the stellar distribution and T^ 6 val- 
ues provide a good description for the stellar components 
of DG1 and DG2. 

From these, we then compute the corresponding stel- 
lar rotation velocities assuming a vertical sech 2 (z) scale 
height distribution of stars. We calculate the vertical 
scale height zq using a ratio of h/zo=2.5 where h is the 
radial scale length of disk, derived from the 3.6/im sur- 
face brightness profile. The derived scale heights z of 
DG1 and DG2 are 0.41 kpc and 0.28 kpc, respectively. 
These are similar to the mean value (0.32 kpc) of the 
7 THINGS dwarf galaxies as given in Oh et al. (2011). 
The derived stellar rotation velocities of DG1 and DG2 
are shown in panels (d) of Fig. 3, and they agree well 
with the true ones indicated by dashed lines. 

3.3. The disk-halo decomposition 

We separate the dark matter components of DG1 and 
DG2 by subtracting the mass models of baryons derived 
in Section 3.2 from their total rotation curves. We then 
fit two halo models, the NFW and pseudo-isothermal 
(ISO) halo models, to these kinematic residuals in order 
to examine the dark matter distribution in a quantita- 
tive way. The NFW and pseudo-isothermal halo models 
represent cusp-like and constant density (core) matter 
distributions at the centers of galaxies, respectively. 

The NFW halo model (Navarro et al. 1996, 1997) is 
given as, 



Vnfw(R) — V2 



do 1 



'ln(l + cx) — cx/(l + cx) 
x[\n(l + c) - c /(l + c)] ' 



(1) 



where c is the parameter quantifying the degree of con- 
centration of the dark matter halo. V200 is the rotation 
velocity at radius -R200 where the density contrast with 
the critical density of the Universe exceeds 200 and x is 
defined as R/R^oo- 

Likewise, the rotation velocity based on the pseudo- 
isothermal halo model is as follows, 



Viso(B) 



\ 



A-KGp R 2 c 



1 — atan [ — 
R \Rc 



(2) 



where po and Rc are the core-density and core-radius of 
a halo, respectively. 

By comparing the fit qualities of these two halo models 
to the kinematic residuals for the dark matter compo- 
nent, we examine which halo model is preferred to de- 
scribe the derived dark matter distributions of DG1 and 
DG2. 

As shown in the upper panels of Fig. 4, compared to 
the CDM NFW halo model, despite its feasible fits, the 
pseudo-isothermal halo model gives a better description 
for the derived dark matter distributions of both DG1 
and DG2 in terms of the fit quality (i.e., reduced % 2 
values). This is also confirmed by the velocity residuals 
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between the dark matter rotation curves and the best 
fitted halo models as shown in the lower panels of Fig. 4. 
The best fitted NFW halo models are too steep to match 
the inner regions of the dark matter rotation curves of 
DG1 and DG2. 

4. THE SIMULATIONS VS. THINGS 

In this section, we compare the derived dark mat- 
ter distribution of DG1 and DG2 with that of the 7 
THINGS dwarf galaxies. THINGS is a high spectral (< 
5.2 kms^ 1 ) and angular (~6") resolution Hi survey for 
34 nearby galaxies undertaken using the NRAO 2 Very 
Large Array (VLA) (Walter et al. 2008). THINGS is 
complemented with other data, such as from the Spitzer 
SINGS survey and ancillary optical B, V and R images 
taken with the KPNO 2.1m telescope (Kennicutt et al. 
2003). These high-quality multi- wavelength data signifi- 
cantly reduce observational uncertainties and thus enable 
us to derive more reliable mass models of the galaxies. 

Oh et al. (2011) performed the dark matter mass mod- 
eling of 7 dwarf galaxies in exactly the same way as for 
DG1 and DG2 as described in Section 3. Basic proper- 
ties of the galaxies are listed in Table 1. In particular, 
the selected 7 THINGS dwarf galaxies have similar obser- 
vational properties as DG1 and DG2, such as resolution 
(—150 pc at the distance of —4 Mpc), maximum rotation 
velocity (< 80 kms -1 ), dynamical mass (~10 9 M ) and 
scale height (-0.3 kpc). In addition the THINGS dwarf 
galaxies are not satellites and they have most likely only 
weakly interacted with larger systems, as are the cases 
of DG1 and DG2. Therefore, this allows us to make a 
direct comparison between the simulations and observa- 
tions, and examine if the simulated dwarf galaxies are 
realistic compared to the dwarf galaxies in the local uni- 
verse. 

4.1. The relation between M s t a r and A/haio 

By comparing the stellar masses of DG1 and DG2 to 
their halo masses, we examine whether the star formation 
efficiency of DG1 and DG2 is comparable with that of 
real galaxies which have similar dynamical masses. This 
is important as the number of stars which form in small 
halos can put strong constraints on baryonic feedback 
and its effects. 

It has been found that galaxies of smaller or larger 
halo masses with respect to the Milky Way appear 
to have inefficient star formation (Navarro & Steinmetz 
2000; Governato et al. 2007; Li & White 2009). The 
low star formation efficiency in high mass galaxies is of- 
ten attributed to AGN feedback (Ciotti & Ostriker 2001; 
Benson & Madau 2003), while the effects of supcrnovae 
are often invoked to explain the low star formation effi- 
ciency in low mass galaxies (Larson 1974; White & Rees 
1978). Recently, Guo et al. (2010) derived the rela- 
tion between halo mass and stellar mass from abun- 
dance matching by combining the stellar mass func- 
tion from the Sloan Digital Sky Survey data release 7 
(SDSS/DR7; Li & White 2009) with the halo/subhalo 
mass function from N-body ACDM simulations (Mille- 
nium and Millcnium-II simulations; Springel et al. 2005; 

2 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



9 




10 9 10 10 10" 10' 2 

M h , (M_) 

halo ©' 

FlG. 5. — The A-f 3 tar -A/ halo relation of DG1, DG2 and the 
THINGS dwarf galaxies (Oh et al. 2011) as well as galaxies from 
the literature (van den Bosch et al. 2001; Stark et al. 2009). The 
Local Group galaxies, LMC (Mastropietro et al. 2005; Guo et al. 
2010), SMC (Stanimirovic et al. 2004; Guo et al. 2010), NGC 6822 
(Valenzuela et al. 2007) and NGC 3109 (Valenzuela et al. 2007) are 
also overplotted. The crrorbars for Mhaio and M s t a r of the LMC and 
SMC arc computed based on the different estimates given by Guo et al. 
(2010) and the other two papers (i.e., Mastropietro et al. 2005 and 
Stanimirovic et al. 2004). The crrorbars for NGC 6822 and NGC 3109 
come from the different mass models that Valenzuela et al. (2007) con- 
sidered for these galaxies to reproduce their kincmatical and photomct- 
rical properties assuming that they are hosted in CDM halos. The solid 
curve is from abundance matching by combining the stellar mass func- 
tion from the SDSS/DR7 with the halo/subhalo mass function from 
the Millenium and Millcnium-II simulations (Guo et al. 2010). As de- 
scribed in Guo et al. (2010), the relation below the stellar mass 10 8 3 
was extrapolated assuming constant slope as indicated by the dashed 
line. See Section 4.1 for more details. 



Boylan-Kolchin et al. 2009; see also De Lucia & Blaizot 
2007 based on semi-analytic models). Similar meth- 
ods and results are found in Conroy & Wechsler (2009), 
Moster et al. (2010) and Trujillo- Gomez et al. (2010). 

Guo et al. (2010) have shown that the star formation 
efficiency in most recent hydrodynamical, cosmological 
galaxy formation simulations is higher than that pre- 
dicted from the M sta r— A/halo relation. In particular, as 
shown in Fig. 5, the stellar masses of DG1 and DG2 are 
about an order of magnitude larger than those inferred 
from the relation. 

However, since the SDSS/DR7 only covers the stel- 
lar mass range from 10 s 3 to 10 118 , the M sta r— -Mhaio re- 
lation outside this range was extrapolated assuming a 
constant slope as indicated by the dashed line as shown 
in Fig. 5 (see Guo et al. 2010). The stellar masses of 
DG1 and DG2 fall within the lower extrapolated region. 
Here, we make a more direct comparison between the 
stellar mass to halo mass ratio between simulated and ob- 
served galaxies. To this end, we compare DG1 and DG2 
with nearby low-mass galaxies from van den Bosch et al. 
(2001), Stark et al. (2009) and the local group as well as 
the 7 dwarfs from THINGS Oh et al. 2011). 

For the halo mass of the THINGS dwarf galaxies, we 
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FlG. 6. — Left: The rotation curve shape of DGl and DG2 as well as the 7 THINGS dwarf galaxies. The dark matter rotation curves (corrected for 
baryons as shown in Fig. 4) arc scaled with respect to the rotation velocity Vb.3 at R0.3 where the logarithmic slope of the curve is dlogV / diogR — 0.3 
(Hayashi & Navarro 2006). The small dots indicate the NFW model rotation curves with V200 ranging from 10 to 90 kms -1 . Sec text for further 
details. The best fitted pseudo-isothermal halo models (denoted as ISO) arc also ovcrplotted. See Section 4.2 for more details. Right: The scaled 
dark matter density profiles of DGl and DG2 as well as the 7 THINGS dwarf galaxies. The profiles arc derived using the scaled dark matter rotation 
curves in the left panel. The small dots represent the NFW models (a~ — 1.0) with V200 ranging from 10 to 90 kms -1 . The dashed lines indicate 
the best fitted pseudo-isothermal halo models (a~0.0). Sec Section 4.3 for more details. 

estimate Mono as follows, 



M 200 [M ]=2OOx 



100 x 



8ttG 

a, 

G 



V: 



200 x 3 



/ ' ZiUU \ 



•1QH, 







;3.29 x 10 5 x 



200: 



(3) 

is the Hubble constant (70.6kms _1 Mpc -1 ; 
2010), G is the gravitational constant (4.3 x 
; x km 2 s- 2 l and V 2 



where Hq 
Suyu ct al 

10~ 3 pc Mg 1 km^ s"^) and V200 in kms -1 is the rotation 
velocity at radius -R200 as given in Eq. 1. However, the 
NFW halo model fails to fit the dark matter rotation 
curves of the THINGS dwarf galaxies, giving negative 
(or close to zero) c values (Oh et al. 2011). To circum- 
vent the unphysical fits, we instead fit the NFW model to 
the rotation curves with only V200 as a free parameter af- 
ter fixing c to 5 which is lower than typical values (e.g., 
8-9; McGaugh et al. 2003) predicted from ACDM cos- 
mology. The fitted V200 values of some galaxies are larger 
than their measured maximum rotation velocities. This 
is because the rotation curves are still rising at the last 
measured points. Moreover, as a larger c value induces 
a smaller V200 an d hence lower halo mass, our choice 
of a low c will provide a robust upper limit for our de- 
rived halo mass, as indicated by the arrows in Fig. 5. As 
shown in Fig. 5, despite the uncertainties remaining in 
these estimates, the stellar masses of DGl and DG2 at 
their given halo masses are consistent with those of real 
galaxies. Both the real galaxies and the simulations de- 
viate from the extrapolated line from the M st ar — -Mhalo 
relation in Guo et al. (2010) at low halo masses. How- 
ever, as discussed in Trujillo- Gomez et al. (2010), there 
still remain uncertainties for dwarfs in the sense that the 
observational data suffer from small number statistics 



and the results of abundance matching are incomplete in 
the low-luminosity tail of the luminosity function. 

4.2. The rotation curve shape 

The rotation curve reflects the total potential (dark 
matter + baryons) of the galaxy and thus it is directly 
related to the radial matter distribution in the galaxy 
(and vice versa). Consequently, the cusp-like dark mat- 
ter distributions in the CDM halos impose a unique shape 
on the rotation curves, which steeply rise at the inner 
regions. Therefore, a relative comparison of galaxy ro- 
tation curves between the simulations and observations 
can serve as a useful constraint for testing the ACDM 
simulations. 

In this context, we compare the rotation curves of DGl 
and DG2 with those of the THINGS dwarf galaxies. In 
order to accentuate their inner shapes, we scale the ro- 
tation curves of both the simulations and the THINGS 
dwarf galaxies with respect to the velocity V0.3 at the 
radius -R0.3 where the logarithmic slope of the curve is 
dlogV/dlogR = 0.3 (Hayashi & Navarro 2006). At the 
scaling radius i?o.3i the rotation curves of both simula- 
tions and the observations are well resolved, which allows 
any differences between them to show up. 

The scaled rotation curves, with the kinematic contri- 
bution of baryons subtracted, are shown in the left panel 
of Fig. 6. We overplot the scaled rotation curves of NFW 
CDM halos (dark-matter-only) with different maximum 
rotation velocities ranging from 10 to 350 kms -1 . We 
choose c values of —9 and —8 for dwarf and disk galaxies 
respectively, which in turn provide V200 values ranging 
from -10 to -90 kms" 1 , and -100 to -350 kms" 1 , 
respectively. Considering that the rotation velocities of 
DGl, DG2 and the THINGS dwarf galaxies at the out- 
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FlG. 7. — The dark matter density profiles of DG1 (left) and DG2 (right). The circles represent the dark matter density profiles derived from 
the dark matter rotation curves shown in Fig. 4. The short dashed and dash-dotted lines indicate the mass density profiles derived using the best 
fitted NFW and pseudo-isothermal halo models in Fig. 4, respectively. The open and filled rectangles indicate the density profiles of DG1LT and 
DG1DM, respectively. Sec Section 2 for more details. The inner slope of the profile is measured by a least squares fit (dotted lines) to the data 
points less than 1.0 kpc as indicated by gray dots. In the right panel the long dashed line indicates a least squares fit, excluding the innermost 
point. The measured inner slope a is shown in the panel. The true dark matter density profiles in the simulations are also ovcrplottcd as indicated 
by open stars. The vertical gray dotted lines indicate the force resolution (86 pc) of the simulations. 



ermost measured radii are all less than ^80 kms -1 , the 
CDM rotation curve with V200~90 kms -1 (i.e., small 
dots in Fig. 6) can be regarded as a lower limit to the 
maximum rotation velocities of the galaxies. 

As shown in the left panel of Fig. 6, the scaled rota- 
tion curves of the THINGS dwarf galaxies are similar to 
each other, showing a linear increase in the inner regions. 
The inner shapes of the rotation curves are better de- 
scribed by pseudo-isothermal halo models (dashed lines) 
than the NFW models. This implies that the THINGS 
dwarf galaxies have core-like dark matter distributions 
at their centers (see Oh et al. 2011 for more discussion). 
Similarly, the scaled rotation curves of DG1 and DG2 
are consistent with those of the THINGS dwarf galaxies. 
They significantly deviate from the CDM rotation curves 
at the inner regions and, like the THINGS dwarf galax- 
ies, they increase too slowly to match the steep rotation 
curves of the CDM halos. 



4.3. The dark matter density profile 

A more direct way to examine the central matter distri- 
butions in galaxies is to convert the galaxy rotation curve 
to the mass density profile. In particular, the measure- 
ment of the inner slope of the profile provides a strin- 
gent observational constraint on the "cusp/core" prob- 
lem. With an assumption of a spherical mass distribu- 
tion for the galaxy halo, the galaxy rotation curve V(R) 
can be converted to the mass density profile p(R) by the 
following formula (see de Blok et al. 2001 and Oh et al. 



2008, 2011 for more details) 
1 



P {R) 



4irG 



VdV 
RdR 




(4) 



where V is the rotation velocity observed at radius R, 
and G is the gravitational constant. Here we do not de- 
contract the halos since in these galaxies adiabatic con- 
traction does not occur and rather expansion happens as 
shown in Govcrnato et al. (2010) (see also Dutton et al. 
2007). 

Using Eq. 4, we derive the dark matter density profiles 
of the THINGS dwarf galaxies, DG1 and DG2 as well as 
the CDM halos whose rotation curves are shown in the 
left panel of Fig. 6. In addition, we also derive the corre- 
sponding mass density profiles of the best fitted pseudo- 
isothermal halo models to the THINGS dwarf galaxies. 
As shown in the right panel of Fig. 6, despite the scatter, 
both DG1 and DG2 have shallower mass density profiles 
than dark matter only simulations. Instead, they are 
more consistent with the THINGS dwarf galaxies show- 
ing near-constant density dark matter distributions at 
the centers. 

In Fig. 7, we compare the derived dark matter density 
profiles of DG1 and DG2 with their true full 3-D dark 
matter density distribution. The inner decrease in the 
actual dark matter density profiles of Fig. 7 is due to 
the shape of the potential in the region below the force 
resolution (86 pc). As shown in Fig. 7, for DG1, the ob- 
servationally derived dark matter density profile robustly 
traces the true values but that for DG2 it is found to be 
on average a factor of three lower than its true value at 
the central regions. This is mainly due to the lower gas 
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rotation velocity of DG2 as shown in panel (f ) of Fig. 3, 
resulting in smaller velocity gradients dV/dR in Eq. 4 
and thus smaller densities. However, considering the un- 
certainties in deriving the profile, the recovered profile is 
acceptable to examine the central dark matter distribu- 
tion. 

We determine the inner density slopes a assuming a 
power law (p ~ r a ) and find them to be a=— 0.31 ± 0.07 
for DG1 and a=-0.49 ± 0.06 for DG2, respectively. If 
we re- measure the slope of DG2, excluding the inner- 
most point which has a large errorbar, the slope is flatter 
(a=— 0.27 ±0.05) as indicated by the long dashed line in 
the right panel of Fig. 7. These slopes deviate from the 
steep slope of ~— 1.0 from dark-matter-only cosmological 
simulations. The profiles of both DG1 and DG2 deviate 
from NFW models beyond about 10 times the force reso- 
lution. This tells us that the baryonic feedback processes 
in dwarf galaxies can affect the dark matter distribution 
in such a way that the central cusps predicted from dark- 
matter-only simulations arc flattened, resulting in dark 
matter halos characterised by a core, as found in normal 
dwarf galaxies in the local universe. 

5. CONCLUSIONS 

In this paper, we have compared the dark matter distri- 
bution of the dwarf galaxies from a novel set of SPH+N- 
body simulations by Governato et al. (2010) with that 
of 7 THINGS dwarf galaxies to address the "cusp/core" 
problem in ACDM. The simulations were performed in 
a fully cosmological context, and include the effect of 
baryonic feedback processes, particularly strong gas out- 
flows driven by SNe. Both the simulated and the ob- 
served dwarf galaxies have similar kinematic properties, 
and have been analyzed in a homogeneous and consistent 
manner as described in Oh et al. (2011). The techniques 
used in deriving dark matter density profiles were found 
to provide accurate results when compared with the true 
underlying profiles, supporting the veracity of the tech- 
niques employed by observers. Therefore, this provides a 
quantitative comparison between the simulations and the 
observations, and allows us to examine how the baryonic 
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feedback processes affect the dark matter distribution at 
the centers of dwarf galaxies. 

From this, we test the general predictions from ACDM 
simulations: (1) the steep rotation curve inherent in the 
central cusp, and (2) the steep inner slope of ~ — 1.0 of the 
dark matter density profiles. We find that the dark mat- 
ter rotation curves of the newly simulated dwarf galax- 
ies rise less steeply at the centers than those from dark- 
matter-only simulations. Instead, they are more consis- 
tent with those of the THINGS dwarf galaxies. In ad- 
dition, the mean value of the inner density slopes a of 
the simulated dwarf galaxies is ~— 0.4 ± 0.1. Compared 
to the steep slope of ~— 1.0 predicted from the previous 
dark-matter-only simulations (including our simulations 
run with DM only), these flat slopes are in better agree- 
ment with a=-0.29±0.07 found in the 7 THINGS dwarf 
galaxies analysed by Oh et al. (2011). 

In conclusion, the results described in this paper con- 
firm that energy transfer and subsequent gas removal in 
a clumpy ISM have the net effect of causing the central 
DM distribution to expand, while at the same time lim- 
iting the amount of baryons at the galaxy center. By the 
present time the DM central profile in galaxies DG1 and 
DG2 is well approximated by a power law with slope a 
of^— 0.4±0.1. These values of a are significantly flatter 
than in the collisionless control run and are in agreement 
with those of observed shallow DM profiles in nearby 
dwarf galaxies. 
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